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This  paper  presents  a  hierarchical  predictive  control  strategy  to  optimize  both  power  utilization  and 
oxygen  control  simultaneously  for  a  hybrid  proton  exchange  membrane  fuel  cell/ultracapacitor  system. 
The  control  employs  fuzzy  clustering-based  modeling,  constrained  model  predictive  control,  and  adaptive 
switching  among  multiple  models.  The  strategy  has  three  major  advantages.  First,  by  employing  multiple 
piecewise  linear  models  of  the  nonlinear  system,  we  are  able  to  use  linear  models  in  the  model  predictive 
control,  which  significantly  simplifies  implementation  and  can  handle  multiple  constraints.  Second,  the 
control  algorithm  is  able  to  perform  global  optimization  for  both  the  power  allocation  and  oxygen  control. 
As  a  result,  we  can  achieve  the  optimization  from  the  entire  system  viewpoint,  and  a  good  tradeoff 
between  transient  performance  of  the  fuel  cell  and  the  ultracapacitor  can  be  obtained.  Third,  models  of 
the  hybrid  system  are  identified  using  real-world  data  from  the  hybrid  fuel  cell  system,  and  models  are 
updated  online.  Therefore,  the  modeling  mismatch  is  minimized  and  high  control  accuracy  is  achieved. 
Study  results  demonstrate  that  the  control  strategy  is  able  to  appropriately  split  power  between  fuel  cell 
and  ultracapacitor,  avoid  oxygen  starvation,  and  so  enhance  the  transient  performance  and  extend  the 
operating  life  of  the  hybrid  system. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells  are  clean  and  highly 
efficient  electrochemical  devices  that  convert  hydrogen  directly 
into  electricity.  They  have  been  considered  as  alternative  power 
sources  for  vehicles,  mobile  robots,  backup  power  sources,  etc.  In 
practical  applications,  PEM  fuel  cells  are  usually  arranged  with  aux¬ 
iliary  power  sources,  such  as  batteries  and  ultracapacitors,  to  form 
hybrid  systems. 

The  electric  loads  supplied  by  a  hybrid  fuel  cell  system  may 
frequently  fluctuate.  Abrupt  changes  in  power  may  cause  oxygen 
starvation  in  the  fuel  cell,  may  overcharge  or  overdischarge  the 
ultracapacitor,  and  may  reduce  the  working  life  of  the  system  in 
a  long  term  [1,2].  Therefore,  sophisticated  power  management  and 
oxygen  control  are  necessary. 

Many  studies  have  been  carried  out  on  power  management. 
Jiang  et  al.  [3]  present  an  adaptive  control  strategy  that  adjusts  the 
output  current  set  point  of  the  fuel  cell.  Ferreira  et  al.  [4]  studied 
a  fuzzy  logic  supervisory-based  power  management  strategy  for  a 
fuel  cell/ultracapacitor/battery  combined  electric  vehicle.  Guezen- 
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nec  et  al.  [5]  and  Rodatz  et  al.  [6]  designed  an  optimal  control 
strategy  to  minimize  the  hydrogen  consumption  in  a  hybrid  fuel  cell 
system.  Zhang  et  al.  [7]  proposed  a  wavelet-transform  algorithm 
to  identify  and  allocate  power  demands  with  different  frequency 
contents  to  corresponding  sources  to  achieve  an  optimal  power 
management  control  algorithm. 

In  the  aspect  of  oxygen  control,  Pukrushpan  et  al.  [8,9  ]  developed 
a  mechanistic  model  suitable  for  the  study  of  controls  in  fuel  cell 
systems.  Vahidi  et  al.  [1,10]  used  a  theoretical  model-based  predic¬ 
tive  control  approach  to  manage  current  and  oxygen  so  as  to  avoid 
oxygen  starvation  and  surge  or  choking  of  the  air  compressor  in 
a  fuel  cell/ultracapacitor  hybrid  system.  To  enforce  constraints  on 
the  oxygen  supply  and  protect  the  fuel  cells  from  oxygen  starvation, 
Sun  and  Kolmanovsky  [11]  use  a  robust  load  governor  to  regulate 
the  current  drawn  from  the  fuel  cell.  These  strategies  have  all  proved 
effective  in  avoiding  oxygen  starvation. 

Model  predictive  control  (MPC)  works  in  a  centralized  manner 
for  constrained  control  problems  through  a  multivariable  mini¬ 
mization  [12].  A  key  advantage  of  MPC  over  other  control  schemes 
is  the  ability  to  deal  with  constraints  in  a  systematic  and  straight¬ 
forward  manner  [13].  Multiple  model  adaptive  control  has  been 
widely  utilized  in  improving  the  transient  response  of  systems  with 
boundary  condition  changes  [14]. 

We  have  noticed  that  power  management  and  oxygen  control 
are  typically  investigated  separately.  That  is  to  say,  power  distribu- 
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tion  command  is  determined  first  according  to  load  requirement 
but  without  interacting  with  real-time  fuel  cell  control,  and  then 
fuel  cell  system  control  (in  which  including  oxygen  supply  con¬ 
trol)  is  determined  to  respond  to  the  power  command  signal  while 
not  to  the  load  requirement  directly.  As  a  result,  the  hybrid  system 
control  is  not  global  optimized.  Our  purpose  in  this  paper,  then,  is 
to  describe  a  centralized  multiple  models-based  predictive  control 
scheme  that  addresses  both  power  distribution  and  oxygen  control, 
and  that  can  provide  a  systematic  and  globally  optimized  solution 
under  multiple  constraints. 

The  proposed  control  scheme  is  designed  and  implemented  as 
follows.  First,  characteristics  of  the  hybrid  system  over  its  whole 
operating  range  are  identified  and  expressed  as  multiple  linear 
discrete-time  models  by  employing  the  fuzzy  clustering  technol¬ 
ogy.  Each  model  corresponds  to  a  typical  operating  zone  of  the 
hybrid  system,  and  the  models  are  updated  online  to  cater  for 
parameter  variations  of  the  real  system.  Second,  constrained  MPCs 
are  designed  for  each  model.  Finally,  an  upper-layer  adaptive  switch 
is  designed  to  determine  the  most  appropriate  model  and  to  switch 
the  corresponding  MPC  as  needed.  The  control  scheme  is  aimed  to 
enhance  the  performance  of  the  system,  and  to  protect  the  hybrid 
system  not  only  by  avoiding  oxygen  starvation,  but  also  by  trading 
off  transient  demands  between  the  fuel  cell  and  the  ultracapaci¬ 
tor,  according  to  constraints  and  weighting  matrices  of  the  output 
errors. 

The  paper  is  organized  as  follows.  In  Section  2,  structure  and 
description  of  the  predictive  control  using  multiple  models  are 
introduced.  Section  3  describes  dynamic  modeling  of  the  hybrid 
system.  Controllers  are  designed  in  Section  4.  In  Section  5,  we 
implement  and  discuss  experiment  and  simulation  results.  Con¬ 
clusions  are  given  in  Section  6. 

2.  System  structure  and  description 

We  focus  on  control  of  electric  power  and  of  the  oxygen  supply. 
We  assume  that  the  hydrogen  is  supplied  at  constant  and  appropri¬ 
ate  pressure,  humidity,  and  temperature,  and  that  wave  effects  are 
insignificant.  These  assumptions  should  not  undermine  the  validity 
of  our  work  because  pressure,  temperature  and  humidity  dynamics 
are  much  slower  than  the  fuel  cell  power  dynamics  which  we  study 
in  this  paper  [1]. 

2.1  System  description 

The  hybrid  fuel  cell  system  studied  in  this  paper,  as  shown 
in  Fig.  1,  is  designed  for  an  automobile  application.  The  electri¬ 
cal  outputs  of  two  PEM  fuel  cell  stacks  are  directly  connected  in 
parallel  to  the  propulsion  motor  load.  An  ultracapacitor  bank  is 
also  connected  to  the  load  through  a  bi-directional  DC/DC  con¬ 
verter  to  form  a  hybrid  PEM  fuel  cell/ultracapacitor  system.  The 
ultracapacitor  bank  should  supply  peak  power  and  should  be 
recharged  by  the  fuel  cell.  The  hybrid  system  can  be  quite  effi¬ 
cient  because  the  fuel  cell  directly  supplies  the  bulk  of  the  load 
demand  while  the  DC/DC  converter  only  operates  to  meet  transient 
demands. 

The  distribution  of  power  between  the  fuel  cell  and  the  ultra¬ 
capacitor  depends  on  the  duty  ratio  of  the  DC/DC  converter,  while 
the  oxygen  supply  to  the  PEM  fuel  cell  is  regulated  by  the  voltage 
applied  to  the  air  compressor,  Vcm •  Duty  ratio  of  a  DC/DC  converter 
is  defined  as  the  ratio  of  switch  on  time  interval,  T0n,  to  switching 
period  T,  i.e. 

d=T-f.  a) 

There  exist  two  duty  ratios  in  the  bi-directional  DC/DC  converter. 
One  duty  ratio,  dc,  is  for  charging  the  ultracapacitor,  and  the  other, 
dd,  is  for  discharging  the  ultracapacitor. 


Oxygen  excess  ratio  is  defined  as 


^o2  = 


Wo  2, in 
Wo 2, ret  ’ 


(2) 


where  W02in  is  the  flow  of  oxygen  into  the  fuel  cell  and  Wo2jct 
is  the  mass  of  oxygen  reacted  in  the  fuel  cell  and  is  related  to  the 
current  drawn  from  the  fuel  cell. 

The  state  of  charge  is  usually  defined  as  the  ratio  of  energy  stored 
in  the  ultracapacitor  to  the  rated  energy  capacity  of  the  ultracapac¬ 
itor  [15],  i.e. 


50C  = 


Yl 

V2 

v  c,  max 


0) 


where  Vc  and  Vc,max  are  the  instantaneous  voltage  and  maximum 
voltage  of  the  ultracapacitor,  respectively.  Vc,max  is  constant  and 
controlling  50C  is  equivalent  to  controlling  Vc/VCimax.  As  a  result, 
we  define  a  new  state  of  charge,  SOC,  which  is  easier  to  control  [1  ], 
i.e. 


SOC  = 


Vc 

i/c,max 


(4) 


2.2.  Control  structure  and  principle 

The  framework  of  the  multiple  model  predictive  control  is 
presented  in  Fig.  2.  It  has  four  major  blocks,  namely  model  pre¬ 
dictive  controllers,  models,  adaptive  switch,  and  the  controlled 
system. 

There  are  n  linear  models  and  corresponding  MPCs  in  the  sys¬ 
tem.  Namely,  Mode//  (i  =  1,  2,  . . .,  n)  is  the  ith  linear  model  of  the 
hybrid  system  at  a  typical  operating  zone  and  MPQ  is  a  model  pre¬ 
dictive  controller  designed  for  Mode//.  Mode/on//ne  denotes  the  unit 
that  updates  each  model  online  according  to  real-time  data  so  as  to 
cater  for  time  variation  of  the  system. 

The  adaptive  switch,  also  called  an  upper-layer  controller,  is 
a  decision  unit  that  determines  the  most  appropriate  MPC  to 
control  the  hybrid  system  during  each  control  period.  Briefly,  its 
operating  principle  is  as  follows.  Firstly,  control  signals  are  fed 
into  the  hybrid  system  and  into  each  model  simultaneously.  Sec¬ 
ondly,  the  output  of  each  model  is  compared  to  the  actual  output, 
respectively.  Finally,  according  to  the  errors,  the  adaptive  switch 
periodically  chooses  the  model  that  best  matches  the  observed  per¬ 
formance  and  then  switches  on  the  corresponding  MPC.  At  each 
period,  only  one  MPC  is  in  service  and  occupies  computational 
resources. 

The  system  design  consists  of  two  major  steps:  identification  of 
the  hybrid  system  and  design  of  the  control.  System  identification 
consists  of  collecting  input  and  output  data  from  the  hybrid  system 
and  dividing  it  into  several  sets  based  on  fuzzy  clustering.  Then 
models  are  identified  from  each  data  set.  Control  design  entails 
design  of  the  MPCs  and  design  of  the  adaptive  switch.  In  design¬ 
ing  the  MPCs,  we  will  focus  on  handling  multiple  constraints  and 
enhancing  computational  speed.  In  design  of  the  adaptive  switch, 
it  is  important  to  design  performance  evaluation  function  and  the 
switch  mechanism. 

3.  Modeling  of  the  hybrid  system 

The  hybrid  system  is  a  multiple  input  and  multiple  output  non¬ 
linear  system.  This  section  establishes  the  bank  of  linear  models 
that  describe  the  hybrid  system.  Input  and  output  data  are  divided 
into  multiple  sets  through  fuzzy  clustering  firstly  and  then  the 
model  for  each  data  set  is  identified.  In  addition,  the  models  will 
be  updated  online  during  utilization  phase. 
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Fig.  1.  Schematic  of  the  hybrid  PEM  fuel  cell  system. 


3.1.  Fuzzy  clustering 

To  characterize  the  hybrid  system  using  piecewise  linear  mod¬ 
els,  fuzzy  clustering  is  chosen  to  divide  the  input  and  output  data 
into  multiple  sets  corresponding  to  linear  areas.  It  contains  three 
steps:  determine  input  and  output  variables,  collect  input  and  out¬ 
put  data,  and  classify  data  using  fuzzy  clustering  method.  The  input 
and  output  variables  are  first  determined  as  follows. 

Input  variables  are  chosen  as 

U(k)  =  ( Vcm  dd  dcf,  (5) 


Output  variables  are  chosen  as 

y(k)  =  ( Xo2  Wcp  psm  Va  1st  lc  SOC)T,  (6) 

where  Wcp  is  air  flow  of  compressor,  psm  is  pressure  of  air  supply 
manifold,  Vst  is  voltage  of  the  fuel  cell,  Ist  is  current  of  the  fuel  cell, 
and  lc  is  current  of  the  ultracapacitor.  Power  demanded  by  the  load, 
Pd,  is  viewed  as  a  disturbance  to  the  system.  The  hybrid  system  can 
now  be  characterized  by  these  input  and  output  variables. 

Zi  =  (yT(i)  uT(i)  )T,  i  =  0,1,2,  N- 1,  (7) 


Fig.  2.  Framework  of  multiple  model  predictive  control. 
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where  N  denotes  the  total  number  of  sampled  data.  u(i)gR3, 
y(z)  e  R7,  and  z(z)  e  R10. 

After  the  input  and  output  data  of  the  system  are  collected, 
a  fuzzy  clustering  approach  with  a  fuzzy  covariance  matrix  [16] 
is  applied  to  classify  the  data  into  n  sets.  The  kernel  issue  of  the 
method  is  to  minimize  the  cost  function 

N- 1  n  N- 1  /  n  \  n 

jf = |  - i )  +  y^Ao^/i - 

i=0  j=l  i=0  \j=l  /  j=l 

(8) 

where  &>Zj  is  a  membership  function  of  sample  zz  to  the  jth  class,  a 
is  a  smoothing  parameter,  dy(Oj)  denotes  the  distance  from  point  Z/ 
to  the  center  v-}  of  the  jth  class  and 

dij(0jk) )  =  (z;  -  V] )TMj(Zi  -vj),  (9) 

Mj  g  R10xl°  is  a  symmetric  and  positive  definite  weighting  matrix, 
Oj  =  (Vj ,  Mj),  the  superscript  k  represents  the  step  number  of  the 
iterative  computing  cycle,  A,z  and  are  Lagrange  multipliers,  pj  is 
a  constraint  of  \Mj\. 

Given  data  {zz-}  and  an  initial  guess  0^\  the  classification  is  pro¬ 
ceeded  as  the  following  iterative  algorithm. 

For  k=  1,  2,  . . .,  use  the  computing  results  of  the  (/<  -  1  )th  cycle 
to  compute  the  variables  of  the  kth  cycle: 


(i)  Compute  the  distance  from  point  zz  to  the  jth  class  according 
to  Eq.  (9). 

(ii)  Compute  co^ 


!/(«-!)' 


-1 


If  dih  =  0  for  some  h,  set  a)ih  =  1,  con  =  0  V/  /  h 
(iii)  Compute  new  estimates 


VN_16U“Z- 
(k)  _  2^j= 0  MjjZi 

j  v^N-1  , 


Elio  < 


M<k>  = 

J 


PjIPjl 


V"  -  Vj)(Zi  -  Vj)T 


E«  < 


(10) 


(11) 


(12) 


Let  k  =  k  + 1  and  recycle  to  (i)  until  a  specified  convergence 
criterion,  A,  is  satisfied. 


Fig.  3.  Flow  chart  of  fuzzy  clustering  algorithm. 


Then  Vi  e  {0, 1, . . .,  N  -  1 },  zz  can  be  classified  to  jth  class,  namely 

Z,  6  j/}l/  =  argmax(®i()j  ,  (13) 

where  Tj  denotes  jth  class  sample  set. 

A  flow  chart  of  the  algorithm  is  shown  in  Fig.  3. 


is  unknown  and  is  estimated  by  the  recursive  least  squares 
algorithm  as 


Gi(k)  =  Gi(k-l)  + 


a,X,(k)e[(k) 

1  +Xj(k)Xi(ky 


(16) 


3.2.  Identification 

The  order  of  the  system  is  determined  through  the  theoreti¬ 
cal  model.  For  each  data  set,  we  use  the  following  linear  model 
to  describe  the  relationship  between  the  input  and  output: 

yi(k)  =  0jXi(k)  Vie {1, 2, ,  n(,  (14) 

where  subscript  i  denotes  the  zth  model,  (9/  e  R(na+nb+Vx7  is  a  sys¬ 
tem  parameter  matrix  and 

Xi(k)  =  (yJ(k-  1),  yT(k  -  tia),  uT(k),  uT(k  -  nb))T . 


fl  if  ||ei(fe)||  >M 

ai={  .  ,  (17) 

I  0  otherwise 

where  denotes  estimation  of  (9;,  M  represents  the  maximum 
error  that  can  be  accepted,  and  ez(/<)  is  the  output  error  between 
the  zth  model  and  the  practical  system,  i.e. 

ei(k)  =  y(k)  -  0jXi(k).  (18) 

The  identification  results  are  then  expressed  as  discrete  state 
space  equation,  i.e. 


(15) 


x(fc+  1)  =  AiX(k)  +  BjU(k)  +  BdiPd(k), 


(19) 
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y(k)  =  Cjx{k)  +  Diu{k),  (20) 

where  x  e  Rlria  and  Pd  denote  state  variables  and  power  demand, 
respectively.  At  e R7rlax7na ,  BjeR7nax3,  BdieR7rla,  C;GR7nax7, 
Df  e  R7x3  are  parameter  matrices  determined  by  (9Z. 

4.  Controller  design 

4.1.  Model  predictive  controller  design 

Fundamentally,  a  MPC  should  predict  the  output  trajectory  of 
a  process  and  compute  a  series  of  control  actions,  subject  to  con¬ 
straints,  that  will  minimize  the  difference  between  the  predicted 
trajectory  and  some  desired  trajectory.  A  key  advantage  of  the  MPC 
over  other  control  schemes  is  its  ability  to  deal  with  constraints  in 
a  systematic  and  straightforward  manner. 

Generally,  the  MPC  control  process  consists  of  future  states  or 
output  prediction,  scrolling  optimization,  and  control  signal  imple¬ 
mentation.  Now  we  use  the  identified  models  to  design  MPC.  An 
important  step  is  to  set  an  objective  function,  which  is  defined  as 

Nu 

min  J  =  V  ((y(k+j)-yAk+j))TQ.(y(k+j)-yr(k+j)) 
u(k),...,u(k+JVu- 1)  V 

j= 1 

+  uT(k  +j)RuT{k  +j))  ,  (21) 

where  Nu  is  predictive  horizon,  y{k+j)  is  the  estimated  out¬ 
put  of  the  system  at  instant  k+j  through  models  based  on 
information  available  at  instant  k.  yT(k+j)  is  the  desired  output 
at  instant  k+j ,  and  QgR7x7,  RgR3x3  are  weighting  matri¬ 
ces  on  output  errors  and  control,  respectively.  We  choose  the 
control  horizon  to  be  equal  to  the  prediction  horizon,  and 
define  Q  =  diag  ( 0xO2  Qwcp  Qpsm  Qyst  Qist  Qic  Qsoc )  and 
R  =  diag  ( Rvcm  Rdd  Rdc ),  where  ,  Qwcp,  Qpsm,  Qvst,  Qist,  Qjc , 
and  Qsoc  are  penalties  on  errors  in  Aq2,  Wcp,  psm,  Vsu  ht ,  Ic,  and 
SOC,  respectively.  Rycm,  Rdd,  and  Rdc  are  penalties  on  VCm,  dd  and 
dd,  respectively.  We  set  QWcp  =  0,  Qpsm  =  0,  QVst  =  0,  Qjst  =  0  and 
Qjc  =  0  so  that  the  corresponding  outputs  are  not  penalized  in  the 
performance  index  and  are  only  used  for  checking  the  constraints. 

Remark.  The  penalty  on  error  in  SOC,  Qsoc.  is  set  to  zero  when 
power  demand  changes  rapidly.  In  the  case  of  increasing  of  power 
demand,  it  is  acceptable  that  the  value  of  SOC  is  less  than  its 
expected  value  so  that  the  ultracapacitor  can  supply  more  for  power 
peaks.  However,  if  Qsoc  =/=  0,  according  to  the  objective  function, 
SOC  will  converge  to  its  desired  value  and  then  the  ultracapaci¬ 
tor  will  supply  less.  Similarly,  in  the  case  of  decreasing  of  power 
demand,  if  Qsoc  0,  less  energy  will  be  charged  into  the  ultraca¬ 
pacitor.  When  power  demand  is  stable,  Qsoc  is  set  to  a  non-zero 
value  so  as  to  adjust  the  value  of  SOC  to  its  desired  value.  SOC  is 
always  also  constrained  by  its  lower  and  upper  limits. 

Substituting  state  Eqs.  (19)  and  (20)  in  (21),  the  equation  is 
abbreviated  as 

min J  =  \UT£2U  +  ( yT(k)L  +  G)U,  (22) 

where  £2  eR3Nux3Nu,  LeR7x3Nu,  GeR1x3Nu  are  constant  matrices 
calculated  through  the  system  model  and  weighting  matrices  Q,  R. 

U  =  (uT(k),  uT(k+  1),  uT(k  +NU  -  1))T.  (23) 

4.1.1.  Constraints 

One  of  the  most  important  advantages  of  MPC  is  its  ability  to 
deal  with  constraints.  In  the  hybrid  fuel  cell/ultracapacitor  system, 
several  constraints  must  be  observed. 


Rapid  variation  on  current  will  harm  the  fuel  cell,  so  it  is  desir¬ 
able  to  constrain  the  rate  of  change  of  fuel  cell  current,  i.e. 

—  AImax  C  ht{k  +  1)  —  lst{k )  <  AImax,  (24) 

where  A Imax  is  the  upper  limit  of  the  rate  of  change,  which  is  deter¬ 
mined  by  characteristics  of  the  fuel  cell. 

The  ultracapacitor  must  not  be  overcharged  or  overdischarged. 
The  state  of  charge  of  the  ultracapacitor  is  an  important  perfor¬ 
mance  index  of  the  hybrid  system,  and  it  should  be  limited  to  some 
desired  range 

SOCmin  —  SOC  <  SOCmax ,  (25) 

where  SOCmin  and  SOCmax  are  the  lower  and  upper  limits  of  SOC, 
respectively.  In  this  paper  the  SOC  is  limited  between  0.45  and  0.95, 
and  the  desired  value  is  0.70. 

Similarly,  the  ultracapacitor  current  should  not  exceed  its  limi¬ 
tations  as 

—Ic,max  —  h  S  Ic,max »  (26) 

where  Ic,max  is  the  upper  limit  of  Jc. 

Starving  the  fuel  cell  of  oxygen  causes  serious  life-limiting  prob¬ 
lems,  while  supplying  it  with  too  much  oxygen  decreases  system 
efficiency.  As  a  result,  the  oxygen  excess  ratio  should  be  limited  as 

^min  —  ^02  —  hmax,  (27) 

where  Xmin  and  Xmax  are  the  lower  and  upper  limits  of  oxygen  excess 
ratio,  respectively.  The  value  of  Aq2  depends  on  oxygen  concentra¬ 
tion  in  the  cathode  and  current  drawn  from  the  fuel  cell.  Generally, 
the  best  value  of  oxygen  excess  ratio  is  2.  In  this  paper  the  range  of 
Xq2  is  constrained  between  1.9  and  2.1. 

The  voltage  of  the  fuel  cell,  also  the  bus  voltage,  is  constrained 
as 

Vst,min  —  Vst  —  ^st,max »  (28) 

where  Vstjnin  and  Vst,max  are  the  lower  and  upper  limits  of  Vst,  respec¬ 
tively. 

The  duty  ratio  of  the  DC/DC  converter  must,  by  definition,  satisfy 
the  following  inequalities: 

0<dd<l.  (29) 

0  <  dc  <  1 .  (30) 

The  voltage  applied  to  the  compressor  motor,  Vcm ,  cannot  be 
larger  than  the  rated  voltage 

0  <VCm<  W  (31) 

Deduced  from  Eqs.  (19)-(21)  and  inequalities  (24)-(31),  the 
control  optimization  is  transformed  to  the  following  constrained 
quadratic  programming  problem 

min I  =  ll/TJ2ir  +  (yT(k)I  +  G)U 

u  2  ^  (32) 

subject  to  Umin  <EU  <  Umax 

where  Umaxe  R8,  and  EeR8x3Nu  are  constant  matrices 

obtained  from  Eqs.  (19)  and  (20)  and  inequalities  (24)-(31). 

4.1.2.  Optimization 

The  constrained  quadratic  programming  problem  is  then  trans¬ 
ferred  to  the  following  simple  optimization  problem  [17] 

KU  +  q  =  F{I<U  +  q  -  0),  (33) 

where  U  g  R8  is  a  variable  to  be  optimized  and 
q  =  -EC2~^{yT(k)L  +  G)T1  K=E£2~^ET,  F  is  a  piecewise  function 
of 

F(U)  =  (F(U1)...F(U8)),  (34) 
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and  for  i=  1, . . .,  8 


F(Ui)  = 


Ui,mim  <  ^i,min 
Ui  5  Ui,min  —  Uj  <  kJ[max 
Ui,maxi  Uijnax  < 


(35) 


where  L/j  min  and  L/j  max  are  the  minimum  and  maximum  values  of 
t/j,  respectively. 

If  D*  is  a  solution  of  Eq.  (33),  then  the  optimal  solution  of  the 
quadratic  programming  problem  (32)  is 

U*  =  &U*  +  (36) 

where  ^  ,  #  =  -1f2_1(yT(/<)L  +  G)T. 

LI  can  be  optimized  through  the  structure  shown  in  Fig.  4. 

MPCs  are  designed  for  each  linear  model,  but  only  one  will  oper¬ 
ate  online  during  each  period.  This  optimization  method  can  also 
be  implemented  through  neural  network  and  it  has  a  faster  conver¬ 
gence  rate  than  related  optimization  methods. 


4.2.  Adaptive  switch  design 

We  describe  next  the  adaptive  switch  that  defines  which  MPC 
operates  online  during  each  control  period.  We  define  a  perfor¬ 
mance  criterion  function  of  the  ith  model  as 

k 

Js,i  =  yefW  +  n  ^2  e-Tr(k-J)e?0’)>  «'  =  1,  •  ••,  n,  (37) 

j=k-l+l 

where  y,  /x,  r  are  constants  and  /  is  an  integer,  e^k)  denotes  the 
error  between  output  of  the  ith  model  and  output  of  the  practical 
system,  i.e. 

ei(k)  =  y,-(fc)  -y{k).  (38) 

The  adaptive  switch  controller  periodically  compares  the  per¬ 
formance  criterion  function  and  switches  on  the  MPC  designed  for 
the  jth  model. 

n 

j  =  argminjs (39) 
i  «'= i 

The  design  procedure  and  operating  principle  of  the  proposed 
modeling  and  control  method  can  be  summarized  as  the  following: 

(I)  Classify  the  real-world  data  according  to  the  flow  chart  in 
Fig.  3. 

(II)  Identify  parameters  <9*  using  Eq.  ( 16)  and  establish  state  space 
models  (19)  and  (20). 

(III)  Choose  a  model  based  on  (39). 

(IV)  Optimize  0  according  to  Fig.  4. 

(V)  Compute  control  variable  U*  using  (36). 

(VI)  Control  the  system. 

(VII)  If  there  is  no  power  demand,  classify  the  new  collected  data 
and  update  models  of  the  system,  else  go  to  (III)  and  start  next 
control  period. 


Table  1 

Parameters  used  in  the  experiment  and  simulation. 


Sym. 

Meaning 

Value 

Tst 

Temperature  of  fuel  cell 

338  K 

Tatm 

Atmospheric  temperature 

288  K 

Ph2 

Partial  pressure  of  hydrogen 

1.5  atm 

n 

Number  of  cells  in  each  stack 

260 

A 

Active  area  of  fuel  cell 

140  cm2 

C 

Capacitance  of  ultracapacitor 

86  F 

Vc,maX 

Maximum  voltage  of  ultracapacitor 

168  V 

5.  Experiment  and  simulation 

The  hybrid  system,  as  shown  in  Fig.  1 ,  is  designed  to  power  a  car. 
The  rated  power  is  50  kW.  The  DC  bus  voltage  is  controlled  between 
170  V  and  250  V.  The  PEM  fuel  cell  is  composed  of  two  PEM  fuel  cell 
stacks,  which  are  connected  in  parallel.  Each  stack  has  260  cells  and 
an  active  area  of  140  cm2. 

The  ultracapacitor  is  86  F  and  the  maximum  voltage  is  168  V. 
The  value  of  capacitance  can  be  realized  by  a  bank  of  56  stan¬ 
dard  ultracapacitors,  each  with  capacitance  of  4800  F  and  a  rated 
voltage  of  3  V,  connected  in  series.  The  maximum  stored  energy 
is  304  Wh  (1/2  x  86  x  1682/3600  x  0.952),  although  only  97  Wh 
(1/2  x  86  x  1682/3600  x  (0.702-0.452))  is  available  between  the 
target  nominal  SOC  and  the  lower  limit.  This  97  Wh  corresponds 
to  an  average  power  at  35  kW  for  10  s  and  that  is  sufficient  to  buffer 
the  fuel  cell  from  acceleration  transients. 

The  experiments  and  simulations  described  in  this  section  relied 
on  Matlab/Simulink  and  were  developed  in  two  steps.  First,  models 
of  the  hybrid  system  were  developed  from  real-world  data  from 
the  system.  Second,  multiple  models-based  MPCs  were  designed 
to  control  the  hybrid  fuel  cell  system  and  comparisons  between 
constrained  and  unconstrained  MPCs  were  made. 

5.1.  Modeling  experiment  and  simulation 

Real-world  input  and  output  data  of  a  50-kW  PEM  fuel  cell  was 
collected.  Assuming  hydrogen  pressure,  humidity,  and  temperature 
of  the  stacks  are  constant.  Related  parameters  used  are  shown  in 
Table  1. 

The  sampled  data  are  equally  divided  into  two  groups.  One  group 
is  used  for  modeling  and  the  other  is  used  for  validating.  In  mod¬ 
eling  phase,  we  first  classify  the  real-world  data  into  5  sets  using 
fuzzy  clustering,  and  then  we  identify  5  models  according  to  each 
data  set  based  on  Eqs.  (16)-(18).  In  each  control  period,  the  most 
suitable  model  is  chosen  from  the  5  models  through  Eq.  (39)  to 
calculate  the  output  of  the  system.  Simulation  result  of  PEM  fuel 
cell  modeling  is  shown  in  Fig.  5.  It  is  shown  that  there  always 
exists  a  model  that  adequately  represents  the  practical  system  and 
the  multiple  models  closely  match  the  practical  fuel  cell.  Multiple 
linear  models  of  other  components  of  the  hybrid  system,  such  as 
DC/DC  converter  and  ultracapacitor  are  established  based  on  their 
mechanistic  models. 


Fig.  4.  Diagram  for  optimization.  F(  )  is  the  piecewise  function  (34),  a  >  0  is  a  scaling  constant. 
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Fig.  5.  Simulation  of  using  the  multiple  models.  5  models  are  used  in  the  modeling,  (a)  Voltage  curves  of  the  PEM  fuel  cell  and  (b)  sampled  current  of  the  PEM  fuel  cell. 


5.2.  Control  simulation 

Control  performances  of  constrained  and  unconstrained  MPCs 
are  studied  and  compared  to  validate  the  proposed  constrained 
MPC.  The  design  of  unconstrained  and  constrained  MPCs  is  similar 
and  is  described  as  the  follows. 

Based  on  multiple  models  established  previously,  MPCs  for  each 
model  and  the  adaptive  switch  are  designed  and  implemented  in 
Matlab/Simulink.  The  output  weighting  matrix  Q.  is  empirically  set 
as  Q=diag(  100,  0,  0,  0,  0,  0,  10),  where  100  is  the  penalty  for  Xq2 
and  10  is  the  penalty  for  SOC.  It  is  noted  here  that  the  penalty  for 
Xo2  is  much  higher  than  that  of  50C.  This  is  because,  in  the  hybrid 
system,  oxygen  excess  ratio  is  the  most  important  index  and  should 
be  controlled  with  high  priority.  The  penalty  for  SOC  is  set  to  zero 
in  case  of  rapid  change  of  power  demand. 

The  input  weighting  matrix,  R ,  which  includes  penalties  for 
air  compressor  voltage  and  duty  ratio  of  the  DC/DC  converter, 
is  set  as  R  =  diag( 0.001,  100,  100).  Since  the  compressor  voltage 
is  usually  over  100  V  and  duty  ratio  is  no  more  than  1,  the  two 
penalties  for  duty  ratios  (100)  are  much  larger  than  that  of  com¬ 
pressor  voltage  (0.001).  In  addition,  control  horizon  is  set  to  10 


Table  2 

Constraints  for  the  constrained  MPC. 


Sym. 

Meaning 

Lower  limit 

Upper  limit 

A  Imax 

Rate  of  change  of  fuel  cell  current 

-20AS-1 

20  As-1 

SOC 

State  of  charge  of  the  ultracapacitor 

0.45 

0.95 

Ic 

Current  of  the  ultracapacitor 

-250A 

250  A 

A-o2 

Oxygen  excess  ratio 

1.9 

2.1 

vst 

Voltage  of  the  fuel  cell 

170  V 

250V 

Vcm 

Voltage  of  the  air  compressor 

0 

200V 

sampling  steps.  The  constraints  of  the  constrained  MPC  are  listed 
in  Table  2. 

A  typical  urban  driving  cycle  in  China  is  used  in  simulation  and 
the  power  profile,  as  shown  in  Fig.  6,  is  considered  as  the  power 
demand. 

The  simulation  results  for  both  the  unconstrained  and  the  con¬ 
strained  MPC  are  shown  in  Fig.  7. 

It  is  noticed  that,  in  Fig.  7(a),  there  exist  significant  perturbations 
in  oxygen  excess  ratio  that  fluctuates  between  1.58  and  2.28  for 
unconstrained  MPC.  As  a  result,  oxygen  starvation  is  unavoidable  in 
this  situation.  While  in  the  case  of  the  constrained  MPC,  the  oxygen 
excess  ratio  is  relatively  stable  and  only  varies  between  1.92  and 
2.07,  which  is  a  much  narrower  range  than  for  the  unconstrained 
MPC.  This  is  because  of  the  constraints  on  oxygen  excess  ratio  and 
the  rate  of  change  of  fuel  cell  current. 

As  shown  in  Fig.  7(b)  and  (c),  current  and  voltage  are  much 
smoother  when  using  the  constrained  MPC  than  when  using  the 
unconstrained  MPC.  The  maximum  rate  of  change  of  the  fuel  cell 
is  20  As-1,  which  is  the  upper  limit  of  the  rate.  The  voltage  curve 
is  inside  the  limited  area  and  stable.  These  phenomena  prove  that 
the  constraint  on  voltage  and  current  change  rate  is  valid.  If  capaci¬ 
tance  of  the  ultracapacitor  increases,  the  upper  limit  of  the  current 
change  rate  can  be  set  lower.  As  a  result,  there  is  a  tradeoff  between 
upper  limit  of  the  change  rate  of  fuel  cell  current  and  capacitance 
of  the  ultracapacitor. 

In  Fig.  7  (d),  the  variation  of  SOC  of  the  ultracapacitor  is  much 
larger  when  using  the  constrained  MPC  than  when  using  the  uncon¬ 
strained  MPC.  This  is  because,  under  the  constrained  MPC,  much 
more  energy  is  squeezed  out  from  the  ultracapacitor  to  supply  the 
peak  load  demands,  and  so  smooth  perturbations  of  fuel  cell  cur¬ 
rent.  The  SOC  satisfies  the  constraints  on  lower  and  upper  limits. 
Therefore,  the  constraints  on  SOC  are  effective.  In  case  of  rapid 
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change  of  power  demand,  penalty  on  50C  in  weighting  matrix  Q. 
is  set  to  zero.  In  this  situation  SOC  is  only  constrained  by  its  lower 
and  upper  limits,  and  the  ultracapacitor  discharges  more  or  charges 
more  energy  to  smooth  the  fuel  cell  transient.  As  a  result,  a  good 
tradeoff  between  the  transient  of  the  fuel  cell  and  ultracapacitor 
can  be  achieved,  according  to  constraints  and  weighting  matrices 
of  the  output  errors. 


The  power  distribution  under  the  constrained  MPC  is  shown 
in  Fig.  8.  It  is  seen  that  the  fuel  cell  power  is  smooth  while  the 
frequent  high  power  demands  are  distributed  to  the  ultracapac¬ 
itor.  The  change  rate  of  the  fuel  cell  current  and  power  are  low. 
Moreover,  main  variables  (e.g.,  X0 2,  SOC,  Ic)  meet  their  constraints. 
Consequently,  the  hybrid  system  is  well  protected,  which  is  helpful 
to  extend  the  lifetime  of  the  system. 


(b)  Current  of  the  fuel  cell 


Fig.  7.  Simulation  results  of  constrained  (red  line)  and  unconstrained  (blue  line)  MPC.  The  weighting  matrices  are  Q=diag(  100,  0,  0,  0,  0,  0, 10),  R  =  diag(0.001, 100, 100).  (a) 
Oxygen  excess  ratio;  (b)  current  of  the  fuel  cell;  (c)  voltage  of  the  fuel  cell  and  (d)  SOC  of  the  ultracapacitor.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  the  article.) 
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Fig.  8.  Power  distribution  of  the  hybrid  system. 


6.  Conclusion 

A  centralized  multiple  models-based  constrained  MPC  for  a 
hybrid  fuel  cell/ultracapacitor  system  was  proposed  and  veri¬ 
fied.  The  proposed  approach,  different  from  other  approaches 
described  in  the  literature,  can  provide  a  globally  optimized  solu¬ 
tion  for  power  distribution  and  oxygen  control,  and  can  protect 
the  hybrid  system  from  oxygen  starvation  by  trading  off  tran¬ 
sient  current  demand  from  the  PEM  fuel  cell  to  the  ultracapacitor. 
It  achieves  this  by  setting  constraints  and  by  weighting  output 
errors. 
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